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Preface 


This  report  was  written  to  provide  the  practicing  engineer  with  a  document 
in  easy  to  understand  terms  and  format  explaining  the  constitutive  model  used 
in  a  nonlinear,  incremental  stmctural  analysis  (NISA).  The  work  was  spon¬ 
sored  under  funds  provided  to  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  by  the  Engineering  Division  of  Headquarters,  U.S.  Army  Corps 
of  Engineers  (HQUSACE),  as  part  of  the  Computer-Aided  Stmctural  Engineer¬ 
ing  (CASE)  Project. 

The  report  was  compiled  and  written  by  Dr.  Kevin  Z.  Tmman,  Washington 
University,  and  Mr.  Barry  D.  Fehl,  Information  Technology  Laboratory  (ITL), 
WES.  The  work  was  managed,  coordinated,  and  monitored  in  the  ITL,  by 
Mr.  Fehl,  Computer-Aided  Engineering  Division  (CAED),  under  the  general 
supervision  of  Mr.  H.  Wayne  Jones,  Chief,  CAED,  and  Dr.  N.  Radhakrishnan, 
Director,  ITL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bmce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors,  Non- 
Si  to  SI  Units  of  Measure¬ 
ment 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

inches 

0.254 

meters 

pounds  per  square  inch  (psi) 

6,894.757 

pascals 

1  Introduction 


In  order  to  perform  a  nonlinear,  incremental  structural  analysis  (NISA)  on  a 
massive  concrete  structure  such  as  a  lock  monolith,  gravity  dam,  or  arch  dam, 
the  structure’s  material  behavior  must  be  accurately  modeled.  This  behavior  is 
modeled  through  the  use  of  a  set  of  constitutive  relationships  that  interact  to 
define  the  material  properties  as  a  function  of  stress,  strain,  temperature,  and 
age  to  be  used  in  the  analysis  to  predict  the  structural  response  to  a  given  load 
history.  The  state-dependent  material  properties  are  developed  for  a  specific 
time  in  the  history  of  the  structure  being  analyzed  and  are  used  to  predict  the 
new  state  of  the  structure  material’s  behavior  at  the  end  of  the  specific  time 
step. 


ANATECH  Model 

When  performing  a  NISA  as  defined  within  ETL  1110-2-365  (31  August 
1994)  requires  the  investigator  (user)  to  provide  a  constitutive  model  for  the 
material.  In  order  to  implement  a  speci^ized  constitutive  model,  a  general 
purpose  finite  element  program  such  as  ABAQUS  with  a  user-defined  material 
subroutine  (called  UMAT  in  ABAQUS)  is  required.  ANATECH  Research 
Corp.  has  developed  a  user-defined  constitutive  model  for  several  different 
applications  including  one  for  mass  concrete  systems  (ANATECH  1992). 
These  constitutive  models  have  been  incorporated  into  a  software  package 
called  ANACAP-U  that  has  been  developed  specifically  for  use  with  the  pro¬ 
gram  ABAQUS. 


Purpose  of  This  Report 

Engineers  without  the  proper  background  or  expertise  would  find  the  devel¬ 
opment  of  a  material  modeling  subroutine  extremely  difficult  if  not  impossible. 
Therefore,  this  report’s  purpose  is  to  discuss  and  define  the  components  of  the 
constitutive  model  to  be  used  within  a  NISA  in  a  direct  and  simplified  manner. 
This  simplified  discussion  should  provide  enough  technical  information  cou¬ 
pled  with  practical  descriptions  to  yield  a  general  understanding  of  modeling 
material  behavior. 
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Currently,  the  constitutive  model  defined  in  ETL  1110-2-365  to  be  used  in 
a  NISA  is  required  to  model  the  concrete’s  modulus,  its  creep  compliance,  its 
autogenous  shrinkage  characteristics,  and  its  cracking  potential.  Each  of  these 
properties  is  dependent  on  several  but  not  all  of  the  following  parameters: 
age,  ambient  conditions,  internal  temperature  distribution,  time,  stress  state, 
strain  state,  loading,  and  concrete  mixture  design  or  concrete  constituents. 
Each  of  these  components  will  be  discussed  in  a  technical  yet  practical 
manner. 


Purpose  of  Model 

The  purpose  of  the  constitutive  model  is  to  accurately  model  the  concrete 
material’s  behavior  during  incremental  constmction  and  its  service  life  in  order 
to  analyze  the  structure's  state  of  stress  and  strain.  For  this  to  occur,  the  strain 
must  be  broken  into  components  that  reflect  the  thermal,  creep,  shrinkage,  and 
mechanical  strains.  The  accuracy  of  these  four  material-related  components 
are  totally  dependent  upon  the  constitutive  model  and  its  relation  to  the  true 
(experimental)  behavior  of  the  actual  material  used  in  construction.  The  strain 
within  a  massive  concrete  structure  is  three  dimensional  and  can  be  broken 
into  four  components  as 

e  =  +  r  (i) 


where 

=  mechanical  strain  (elastic  and  plastic) 

=  thermal-induced  strains 
=  shrinkage  strains 
=  creep-related  strains 

Given  a  state  of  stress  and  strain  at  the  beginning  of  an  increment  and  the 
incremental  total  strain  Ae,  the  constitutive  model  must  compute  the  state  of 
stress  o  and,  for  implicit  calculations,  the  material  tangent  stiffness  matrix  (i.e., 
the  Jacobian  of  the  constitutive  law  dAo/dAe'")  at  the  end  of  the  increment. 

The  stress  and  “tangent”  constitutive  matrix  must  be  compatible  with  the  fail¬ 
ure  and  loading  surfaces  of  the  material  such  as  cracking  in  tension  and  yield¬ 
ing  in  compression.  Since  a  constitutive  law  relates  stresses  to  mechanical 
strains,  the  thermal,  shrinkage,  and  creep  strains  must  be  removed  from  the 
given  total  strains  and  strain  increments  by  the  model  before  defining  the 
stresses.  ANATECH  Research  Corp.  has  developed  such  a  model  for  use  in 
the  ABAQUS  general  purpose  finite  element  code. 
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2  Creep 


After  concrete  is  loaded  it  begins  to  deform  continually  with  time.  This 
deformation  has  two  components:  the  first  is  the  mechanical  deformation 
which  is  immediate  and  the  second  is  time  dependent  and  can  continue  for 
years.  This  second  component  of  deformation  is  referred  to  as  creep.  Creep  is 
commonly  defined  as  the  dimensional  change  or  increase  in  strain  (elongation 
or  shortening)  with  time  due  to  a  sustained  stress.  Creep  can  be  beneficial  or 
detrimental  to  a  structure.  Depending  on  the  structural  configuration,  creep 
can  cause  excessive  deformation  at  later  times,  but  it  can  also  relieve  stresses 
at  locations  of  stress  concentrations  by  redistributing  the  stresses.  For  struc¬ 
tures  with  sustained  thermal  gradients,  creep  can  reduce  the  initial  thermal 
stresses,  but  it  can  also  produce  stress  reversals  once  the  structure  has  cooled 
(Fintel  1974). 

Specific  creep  is  defined  as  the  creep  strain  per  unit  of  sustained  stress. 

The  ultimate  specific  creep  usually  ranges  from  0.2  to  2.0  millionths  per  psi.' 
Very-early-age  deformation  due  to  creep  is  viscoplastic  in  nature  and  predomi¬ 
nantly  unrecoverable,  whereas  deformation  due  to  creep  in  older  concrete  tends 
to  be  viscoelastic  and  predominantly  recoverable.  The  lack  of  creep  strain 
recovery  in  very  young  concrete  is  primarily  due  to  the  rapidly  aging  modulus, 
which  gives  the  material  an  apparent  viscoplastic  behavior.  However,  in 
mature  concrete  aging  is  small  and  creep  strain  recovery  is  almost  total  given 
sufficient  time. 

Controlling  creep  is  a  formidable  task  since  it  is  a  function  of  the  concrete 
constiments,  the  concrete’s  age  at  the  time  of  loading,  the  dimensions  of  the 
structure,  ambient  conditions,  and  curing  procedures.  Constituents  that  affect 
creep  include  the  water-cement  ratio,  aggregate  size,  aggregate  type,  and 
cement  paste  characteristics.  If  only  the  cement  paste  constituents  ate  consid¬ 
ered,  the  higher  the  rate  of  strength  gain,  the  lower  the  creep  deformation  for  a 
given  stress.  Therefore,  an  important  parameter  in  creep  prediction  is  the  ratio 
of  applied  stress  to  the  strength  at  the  time  of  loading.  Creep  is  generally  a 
linear  function  of  the  stress-strength  ratio  up  to  50  percent  of  the  concrete’s 
ultimate  strength  at  which  point  the  creep  becomes  nonlinear.  The  aggregates 


^  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  umts  is  presented  on 
page  vi. 
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affect  creep  by  providing  restraint  with  respect  to  creep  deformation  of  the 
cement  paste.  A  larger  volume  of  aggregate  generally  produces  a  larger  reduc¬ 
tion  in  creep  deformation.  (This  action  also  depends  on  the  physical  properties 
of  the  aggregate  such  as  the  modulus  of  elasticity,  porosity,  grading,  and 
bondability.)  Creep  is  temperature  dependent.  The  creep  strain  can  increase 
for  some  concrete  mixtures  for  temperatures  between  70  and  100  °F  (Fintel 
1974). 

Loading  rate  has  a  significant  effect  on  the  creep.  For  a  step  load,  the  rate 
of  creep  wiU  be  relatively  high  during  the  first  few  days  of  loading,  but  will 
eventually  approach  zero  after  some  finite  amount  of  time.  The  age  of  loading 
has  a  significant  effect  on  the  creep-dependent  deformations.  Early-age-loaded 
concretes  wUl  have  much  larger  ultimate  specific  creep  than  a  specimen  loaded 
at  a  later  age  due  to  the  fact  that  specific  rate  decreases  as  the  strength-to- 
stress  ratio  increases  at  the  time  of  loading.  Figure  1  shows  the  total  (elastic 


Concrete  Age  of  concrete 

plo-cenent  at  tine  of  loading 

with  constant 
stress,  Oo 


Figure  1 .  A  generalized  specific  strain  curve  for  concrete 
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and  creep)  specific  strain  as  a  function  of  time  and  loading  age.  The  elastic 
portion  of  the  curve  (denoted  as  specific  elastic  strain  in  Figure  1)  is  the 
inverse  of  the  aging  modulus.  The  curve  above  the  elastic  response  is  the 
apparent  creep  curve  which  is  greatly  influenced  by  the  age  at  loading.  (See 
discussion  below  for  the  definition  of  apparent  creep  versus  true  creep.) 

Figure  2  gives  typical  experimental  curves  for  a  concrete  cylinder  that 
shows  the  dependence  of  specific  strain  on  the  age  of  the  concrete  at  the  time 
of  loading.  These  experimentally  developed  compliance  curves  (specific 


Figure  2.  Experimentally  derived  creep  strain  curves  for  concrete  cylinders 
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strain)  are  composed  of  two  distinct  components.  The  first  data  point  on  each 
of  the  curves  represents  the  elastic  specific  strain  component.  The  second 
component  is  the  time  curve  which  represents  the  apparent  creep.  The  older 
the  concrete,  the  lower  the  values  of  the  elastic  and  creep  strains. 

Figure  3  is  a  generalized  specific  creep  curve  for  concrete  loaded  with  a 
umaxial  stress  of  Oq  at  age  Tq  and  can  be  used  to  show  the  effects  of  aging  in 
concrete  compliance.  The  first  component  is  the  elastic  specific  strain,  eVoq. 
represented  by  the  vertical  portion  of  the  curve.  This  specific  elastic  strain  can 
be  denoted  as  the  inverse  of  the  elastic  modulus  at  age  Xq,  1/£’(Xo).  At  age  x„ 
the  concrete  has  been  unloaded,  thereby  causing  a  reduction  in  the  specific 
strain  of  l/£(x„).  As  seen  in  Figure  3,  /(Mq)  is  the  measured  data  using  the 
initial  specific  elastic  strain  as  the  reference  datum.  Hence,  /(Mq)  is  called  the 
apparent  creep  compliance.  In  reality  as  the  concrete  ages  to  x„,  the  true  spe¬ 
cific  creep  strain  or  creep  compliance  should  be  Cit-x^,  x„)  which  reflects  the 
changing  modulus.  Therefore,  the  total  specific  strain  can  be  written  in  either 
of  two  forms: 


Figure  3.  Compliance  curve  for  concrete 
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"a 


1 


1 


(2) 


44 


C(t  -  To,T„) 


which  can  be  rearranged  to  give  the  true  creep  compliance  as: 


C(t  -  = 


_4o) 


+  J\ 


('  -  ■'o) 


(3) 


For  nonaging  materials,  ^CXq)  =  and  in  that  case  Equation  3  shows  that 
J  and  C  would  be  equivalent  If  the  concrete  is  mamre  prior  to  loading,  E 
becomes  nearly  constant  and  J  and  C  are  nearly  equivalent  as  shown  in  Fig¬ 
ure  4  for  r  =  X3.  For  early-age  loading,  J  and  C  can  be  quite  different  depend 
ing  on  the  constituents  of  the  concrete.  Quite  often  the  experimentally  derived 
data  for  J  are  used  as  the  creep  compliance  term  causing  some  error  in  the  pre¬ 
dicted  versus  the  experimental  data.  This  error  is  due  to  the  aging  modulus 
and  is  the  difference  between  1/E(Xq)  and  1/E(t„)  as  seen  in  Figure  3  and 
Equation  3.  For  modeling  purposes  separate  identification  of  the  curves  C  and 
J  is  not  necessarily  required.  What  is  required  is  a  model  which  adequately 
represents  the  experimental  behavior  in  loading  and  unloading  requiring  the 
experimentally  derived  specific  strain  curve  and  an  accurate  aging  modulus 
model. 

Creep  compliance  (specific  strain  beyond  the  elastic  specific  strain)  has 
been  modeled  by  ANATECH  in  a  computer  subroutine  as  a  single  equation 
which  is  dependent  upon  time  t,  age  x,  and  temperature  T. 
the  form: 


y(r,x,7)  =  A,(x,T){l  -  e  *  D{x,T)  *  (t 

i=l 


The  equation  is  of 

-  X)  (4) 


where 


A,(x,T) 


£(To) 


£(x) 


(5) 


and 
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Concrete  compliance  curves  at  different  ages  of  loading 


D{%,T) 


E{Xq) 


(6) 


where 


R  =  universal  gas  constant  of  1.98 

•^2’  ^1’  ^2'  P'  Q'  ^  =  mathematically  derived  constants  using  experi¬ 
mental  data 

Xq  =  reference  age  at  which  creep  specimen  is 
loaded 
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These  equations  can  be  used  to  represent  the  experimental  results  by  determin¬ 
ing  the  appropriate  values  for  the  constants  /ij,  A2,  ^2*  Th^  experi¬ 

mentally  derived  constants  required  for  this  equation  ensure  a  consistent  match 
between  the  experimental  data  and  the  computerized  constitutive  model.  The 
constants  are  determined  using  this  procedure: 

Step  1.  Select  an  experimental  creep  curve  for  a  given  age  Tq-  3-day  creep 
is  a  common  age  for  early  creep  experiments. 

Step  2.  Select  three  times  representing  early,  intermediate,  and  long-term 
values  in  order  to  obtain  an  adequate  curve  fit  for  each  region.  Assume  4, 
31,  and  93  days. 

Step  3.  Assume  e  to  be  a  small  value,  0.001  to  0.005,  which  will 

cause  the  first  term  (1  -  e  to  be  dominant  at  early  times.  Solve  for 

rj  (i.e.,  c  0.001  =>  rj  =  6.908). 

Step  4.  Assume  e  to  be  a  small  value,  0.001  to  0.005,  which  will 

cause  the  second  term  to  be  dominant  in  the  intermediate  range.  Solve  for 

r2  (i.e.,  e  o.OOl  =>  r2  =  0.247). 

Step  5.  Select  three  values  for  (r-x),  continue  to  use  x  =  3  days  and  the 
3-day  experimental  creep  curve  and  find  the  appropriate  experimental  creep 
compliance  values.  Construct  these  three  equations: 


y(ri,3)  =  Ai(l  -  e  +  ^2(1  -  e  +  D  ♦  (fj  -  ^ 


J(t2,3)  =  aJi  -  e  +  A2(i  -  e  +  D  *  (^2  -  ^ 


7(r3,3)  =  Ai(l  -  e  +  A2(l  -  e  ^  D  *  (tj  -  ^ 

and  solve  for  Aj,  A2,  and  D.  It  should  be  noted  that  because  of  the  linear 
term  in  Equation  4,  the  creep  rate  at  long  time  approaches  a  constant  value 
instead  of  zero,  which  strictly  is  not  true.  However,  the  error  is  small,  and 
on  the  conservative  side,  for  load  durations  of  interest.  This  situation  can 
be  corrected  by  replacing  the  linear  term  by  another  exponential  term  where 
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/-g  would  be  found  using  Steps  3  and  4  and  ^43  would  be  determined  in 
Step  5.  Steps  6-11  need  only  be  performed  if  T  is  considered  as  a  variable. 

Step  6.  Select  a  reference  temperature  Tq  and  another  temperature  T>Tq 
(Kelvin). 

Step  7.  Determine  the  creep  compliance  values  for  a  fixed  time  t  and  age  x 
(31  and  3  days,  respectively,  are  reasonable  values)  from  experimental 
creep  data  for  temperatures  Tq  and  T. 

Step  8.  Construct  this  equation: 


g-Q/RT 

-Q/RTo 

e 


(10) 


Solve  for  Q. 

Step  9.  If  experimental  creep  data  are  available  for  a  third  temperature  T^, 
calculate  the  specific  creep  strain  y(31,3,7’3)  from  the  fully  defined  equation 
and  compare  with  the  experimental  creep  data  for  the  chosen  time,  age,  and 
ternperamre  (31,  3,  T^,  respectively). 

Step  10.  If  the  error  is  large,  adjust  Q  in  a  manner  which  wiU  reduce  the 
error.  Repeat  until  the  error  is  acceptable. 

Step  11.  Check  for  other  values  of  t,  x,  and  Ty  Adjust  Q  until  there  is 
convergence  or  the  error  is  within  an  acceptable  accuracy. 

After  the  steps  outlined  above  are  performed,  the  resulting  effect  of  the 
constants  computed  must  be  determined  by  performing  analyses  of  a  single 
element  using  the  loading  that  was  applied  in  the  laboratory  for  1-,  3-  and 
14-day  creep  tests.  The  resulting  strains  of  these  analyses  should  be  compared 
directly  with  the  strain  data  obtained  in  the  laboratoiy.  If  the  results  of  all 
three  analyses  do  not  compare  well  with  the  test  data  then  the  constants  must 
be  adjusted,  the  analyses  performed,  and  once  again  the  results  compared  with 
the  test  data.  This  iterative  process  must  continue  until  a  satisfactory  compari¬ 
son  of  all  three  strain  tests  is  achieved.  The  process  used  in  adjusting  the  con¬ 
stants  is  demonstrated  in  Appendix  A. 

Figure  4  shows  different  creep  compliance  curves  for  concrete  loaded  at 
various  ages.  Concrete  loaded  at  age  Xq  has  a  different  compliance  curve  than 
a  concrete  loaded  at  age  Xj.  The  assumption  that  the  reduction  in  the  rate  of 
creep  with  increasing  age  is  proportional  to  the  increase  in  the  modulus  with 
increasing  age  provides  the  reasoning  for  the  use  of  a  ratio  of  the  modulus  at 
different  ages  to  adjust  the  creep  compliance  curve  for  aging.  With  a  reference 
age  of  3  days,  the  ratio  becomes  £'(3)/£(x).  Figure  5  shows  creep  compliance 
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Figure  5.  Concrete  compliance  curves  as  a  function  of  duration  of  loading 


as  a  function  of  load  duration,  t  -  x.  For  equal  durations  of  load,  it  is  shown 
that  early  loading  produces  significantly  more  creep-related  strain  than  a  late- 
age  loading.  The  ratio  of  the  concrete  compliance  curve  values  at  a  specific 
duration  of  loading  is  the  true  adjustment  factor.  From  experimental  data  it 
has  been  determined  that  the  aging  effects  are  best  represented  with  a  factor  of 
[E(3)/E(x)f  where  p  is  between  one  and  two.  As  seen  in  Equations  5  and  6, 
this  factor  is  used  to  account  for  the  aging  effects  in  the  creep  equation.  The 
value  for  p  can  be  experimentally  determined. 
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Concrete  creep  is  assumed  to  follow  a  thermo-itieologicaUy  simple  behavior 
which  allows  the  creep  curve  at  any  temperature  T  to  be  generated  from  a 
reference  creep  curve  developed  for  temperature  Tq  (reference).  This  type  of 
behavior  implies  that  a  creep  curve  at  any  temperature  when  plotted  as  a  func¬ 
tion  of  log  t  can  be  obtained  by  a  simple  shift  of  the  reference  curve  along  the 
log  t  axis.  The  shift  ^(T)  is  called  a  shift  factor  and  the  log  of  (j)(T)  is  called 
the  shift  function.  This  behavior  allows  the  implementation  of  an  extrapola¬ 
tion  scheme  that  can  cover  a  wide  range  of  temperatures  with  only  a  few 
temperature-controlled  creep  tests.  Equation  9  can  be  rewritten  to  adjust  for 
temperature  effects  in  creep  compliance  as: 


.-Q/RT 

J(t,x,T)  -  JiuxjQ)—^— 


(11) 


The  factor  for  adjusting  the  creep  compliance  for  temperature  effects  in  the 
ANACAP-U  subroutine  becomes: 


0(7)  =  c  (12) 


since  the  terms  related  to  Tq  are  already  coupled  in  the  experimental  data  used 
to  generate  the  creep  compliance  constants  Aj,  ^2-  rj,  rj,  and  D.  Once  again 
if  temperature  is  found  to  be  of  httle  or  no  impx)rtance  the  use  of  thermal 
adjustment  factors  is  unnecessary. 
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3  Autogenous  Shrinkage 


Autogenous  shrinkage  is  a  decrease  in  volume  of  a  concrete  specimen  or 
member  due  to  hydration  of  the  cementitious  materials  without  gaining  or 
losing  moisture.  This  type  of  shrinkage  occurs  within  the  central  regions  of 
mass  concrete  systems.  For  mass  concrete  systems  this  component  of  shrink¬ 
age  can  be  significant  compared  to  drying  shrinkage.  Autogenous  shrinkage 
occurs  over  a  much  longer  time  period  than  drying  shrinkage  which  is  a  local¬ 
ized  phenomenon  that  affects  only  a  thin  layer  of  concrete  near  the  concrete 
ambient  air  interface.  Autogenous  shrinkage  increases  with  respect  to  higher 
cement  content  and  increased  temperatures.  Sealed  creep  cylinders  with  no 
external  loading  have  been  successfully  used  as  a  means  of  measuring  autoge¬ 
nous  shrinkage. 

The  autogenous  shrinkage  has  been  modeled  by  ANATECH  in  a  computer 
subroutine  as  a  single  equation  with  four  required  constants  and  two  time- 
dependent  exponential  terms: 


=  Ci(l  -  e  +  €2(1  -  e 


These  constants  are  evaluated  in  a  manner  similar  to  that  for  the  creep 
constants. 


Step  1.  Select  a  late  time  r  (100  to  3(X)  days).  Assume  that  the  term  e 
is  small  (0.001  to  0.005).  Solve  for  S2  (i.e.  e  _  q.oOI  =>  ^2  = 

0.0691). 


-SjClOO) 


Step  2.  Select  an  intermediate  time  t  (20  to  60  days).  Assume  that  the 

term  e  is  small  (0.001  to  0.005).  Solve  for  (i.e.,  e  = 
0.001  =>  =  0.1727). 

Step  3.  Select  two  data  points  from  the  experimental  data  e^(30)  and 
6^(100).  Substitute  these  values  into  Equation  13. 


,(l00).c,(l  -  -q(l  - 
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ej30)  =  Ci(l  -  e  +  cj,!  - 


(15) 


Step  4.  Solve  the  equations  for  the  constants  C,  and  €2- 

Step  5.  Once  the  constants  have  been  determined,  the  equation  using  these 
constants  should  be  compared  with  the  experimental  data  in  order  to 
explore  the  accuracy  of  the  generated  curve.  Adjust  5]  and  ^2  until  a  satis¬ 
factory  curve  is  generated. 
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4  Aging  Modulus 


The  modulus  of  elasticity  for  concrete  is  generally  defined  from  compres¬ 
sive  tests.  It  is  heavily  dependent  upon  the  cement  paste,  aggregate  modulus 
of  elasticity,  and  the  relative  volumes  of  aggregate  to  cement  paste.  The 
dependence  of  modulus  of  elasticity  on  the  cement  paste  is  the  reason  that  the 
modulus  is  also  concrete-age  dependent.  The  cement  paste  modulus  is  depen¬ 
dent  upon  the  amount  of  hydration  that  has  occurred,  increasing  as  the  hydra¬ 
tion  process  continues.  Also  the  increased  bonding  of  the  cement  paste  to  the 
aggregate  improves  the  composite  action  of  the  cement  paste  and  aggregate 
which  is  also  a  function  of  the  hydration  or  curing  of  the  concrete.  Typically, 
the  modulus  of  elasticity  increases  faster  than  the  compressive  strengths  at 
early  ages  (Fintel  1974). 

The  aging  modulus  as  a  function  of  concrete  age  has  been  modeled  in  the 
ANATECH  version  of  UMAT  (ANACAP-U)  by  a  single  equation.  Within  this 
equation  several  constants  must  be  determined  in  order  to  fit  the  curve  to  the 
experimental  aging  modulus  data.  The  equation  is  of  the  form: 

E(x)  =  Eg  +  Ei(l  -  e  +  £2(1  -  e  +  E3T 

where  T  is  the  age  of  the  concrete  given  in  days.  One  way  to  fit  the  equation 
to  the  actual  aging  modulus  versus  age  data  is  to  generate  four  equations  to 
solve  for  the  four  constants  Ej.  Each  term  in  the  equation  is  used  to  best 
represent  a  region  of  the  actual  aging  modulus  data.  Eg  is  used  to  represent 
the  first  day;  E^^  is  used  to  represent  the  data  for  early  values  (5  to  10  days); 
Ej  is  used  to  represent  the  data  for  intermediate  values  (20  to  40  days);  and  E3 
is  used  to  represent  late  data.  Typically  the  constants  would  be  determined  by 
this  procedure: 

Step  1.  Eq  is  the  value  of  the  aging  modulus  at  t  =  1  day.  Any  value  of 
the  aging  modulus  required  prior  to  x  =  1  day  is  determined  using  linear 
interpolation  from  0  to  the  value  for  day  1. 
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Step  2.  Select  an  intermediate  age  such  as  28  days  where  aging  effects  are 
small  or  can  be  neglected.  Assume  the  term  e  is  small  (0.001  to 

0.005).  Solve  for  Cj  (i.e.,  e  =  0.001  =>  =  0.2558). 

Step  3.  Select  an  early  age  such  as  7  days  where  the  intermediate  age 
effects  are  small  or  can  be  neglected.  Assume  the  term  e  is  small 
(0.001  to  0.005).  Solve  for  C2  (i.e.,  =  0.001  =>  C2  =  1.151). 

Step  4.  Determine  the  experimental  values  for  E(7),  £(28),  and  £(90  days 
or  later)  from  the  experimental  data.  The  values  of  7  and  28  days  are  usu¬ 
ally  convenient  since  many  experiments  are  developed  to  collect  data  at 
those  two  concrete  ages. 

Step  5.  Solve  the  corresponding  four  equations  for  £;,  i  =  0...3. 

EO)  =  £0  (17) 


£(7)  =  £0  +  £,(l  -  +  £2(1  -  +£3*7  (IS) 


£(28)  =  £0 

+  £i(] 

1  -e-'®-")  -£,(1  -e->®-’')*£3 

*  28 

(19) 

0 

II 

+  £i(] 

*  90 

(20) 

Step  6.  Once  the  six  constants  have  been  determined.  Equation  16  should 
be  compared  with  the  experimental  data  in  order  to  explore  the  accuracy  of 
the  generated  curve.  Adjust  the  constants  to  obtain  a  reasonably  accurate 
equation.  Again  as  in  the  case  of  continuously  diminishing  creep  rate  at 
long  times,  the  modulus  aging  should  eventually  saturate.  Equation  16, 
however,  implies  continuous  aging  approaching  a  linear  rate  at  long  times. 
If  long-term  loading  is  of  interest.  Equation  1 6  may  have  to  be  refitted  with 
a  fourth  exponential  term  instead  of  a  linear  term. 
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Material  Subroutine 


ABAQUS  allows  the  use  of  a  predefined  constitutive  model  through  the 
ABAQUS  UMAX  subroutine.  ABAQUS  provides  a  warning  that  this  option 
should  only  be  used  by  experts  and  that  any  user-developed  constitutive  model 
should  be  thoroughly  tested  prior  to  its  implementation.  ANATECH ’s  heavy 
use  of  ABAQUS  and  its  UMAX  subroutine  in  the  analysis  of  nuclear  facilities 
provided  the  expertise  needed  in  developing  a  constitutive  model  for  incremen¬ 
tally  constructed  massive  concrete  structures.  ANAXECH’s  development  of 
ANACAP-U  includes  an  early-age  constitutive  model  which  can  be  used  for 
massive  concrete  structures  on  Corps  of  Engineers  projects. 

In  order  to  model  the  complex  behavior  of  creep,  shrinkage,  aging  modulus, 
and  cracking,  a  constitutive  model  subroutine  was  developed  to  be  used  with 
the  ABAQUS  finite  element  code.  Nonlinear  finite  element  solutions  are 
highly  dependent  upon  the  numerical  implementation  of  a  sound  constitutive 
model.  The  constitutive  model’s  purpose  is  to  accurately  define  the  stresses 
and  the  tangent  constitutive  matrix  at  the  end  of  each  time  step.  In  a  finite 
element  program  the  required  equations  would  be  of  the  form: 

=  Ae  -  Ae^  -  -  Ae<^ 


where 

D  =  the  tangent  constitutive  matrix 

6-  =  the  stresses  at  time  i  =  t  ot  t  +  At 

At  =  the  difference  in  the  strains  at  r  -h  At  and  t 

The  increment  of  total  strains  At  in  Equation  21  are  those  strains  that 
result  from  the  solution  of  the  force-displacement  (equilibrium)  equations  and 
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are  returned  by  ABAQUS  to  the  constitutive  subroutine  ANACAP-U.  The 
other  strain  increments  Ag‘  due  to  temperature,  shrinkage,  and  creep  are  cal¬ 
culated  within  the  constitutive  routine.  The  UMAT  subroutine  is  called  by 
ABAQUS  whenever  a  ^MATERIAL  definition  inside  the  user  input  includes  a 
*USER  MATERIAL  option  to  define  the  mechanical  constitutive  behavior  of 
the  material.  The  *USER  MATERIAL  option  forces  every  element  defined  by 
the  *MATERIAL  name  to  be  controlled  by  the  UMAT  subroutine.  The 
*USER  MATERIAL  option  requires  a  number  that  reflects  the  number 
of  user-defined  properties.  The  properties  to  be  supplied  are  defined  in 
the  ANACAP-U  User’s  Manual  (ANATECH  Research  Corp.  1992)  and 
are  presented  in  Table  1.  For  further  information  regarding  the  property 
definitions,  refer  to  the  ANACAP-U  User’s  Manual. 

A  portion  of  the  ABAQUS  input  for  the  Olmsted  Locks  NISA  is 
presented  in  Chapter  6  to  illustrate  the  proper  use  of  the  ABAQUS 
commands  and  the  required  input  for  using  the  UMAT  subroutine  as 
contained  in  ANACAP-U. 
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Table  1 

Material  Properties  Required  for  UMAT  Input 

Prop. 

No. 

Description 

Comments 

1 

Model  Flag 
=1  Elastic  Comp. 

=2  Elastic  Perfectly 
Plastic  Comp. 

=3  Comp.  Strain 

Hard,  and  Soft. 

=4  Elastic  Perfectly 
Plast.  Comp,  with 
Creep,  Aging  and 
Shrinkage 

This  flag  should  be  set  to  4  for  NISA  studies 

2 

Static  crushing 
strength  (f;,') 

Input  f^'  from  concrete  testing  tor  concrete  3  days  old. 

Typical  values  are  600  to  1,000  psi 

3 

Static  tensile 
cracking  strain 

Calculate  using  equation  from  ETL  1 '[10-2-365  and  results 
from  slow  load  beam  test  performed  on  selected  project 
mixture 

4 

Static  elastic  modulus 
(E) 

Input  E  from  concrete  testing  for  concrete  3  days  old. 

Typical  values  are  1,500,000  to  2,500,000  psi 

5 

Poisson's  ratio  (v) 

Input  from  concrete  testing.  Typical  values  are  0.15  to  0.2 

6 

Coef.  of  thermal 
expansion  (a) 

Input  from  concrete  testing.  Typical  values  are  4.0  x  10  ®  to 

6.0  X  10-®  in./in./®F 

7 

Stress  free  temp. 

Use  placing  temperature 

8 

Flag  for  English  or 

SI  Units 

Use  a  1  for  English  units 

9  to 

24 

Properties  9  through  24  are  not  used  in  a  NISA  and  should  all  be  set  to  0.0. 

Information  on  these  properties  can  be  found  in  the  ANACAP-U  User's  Manual. 

25 

Creep  fit  flag 

Use  a  value  of  2 

26 

Age  of  concrete  in 
days 

Typically  this  will  be  0.0 

27 

Shrink  factor 

Multiplier  for  shrinkage  curve,  should  be  based  on  bounds 
being  applied  to  shrinkage  as  specified  in  ETL  1 110-2-365. 

If  15%  increase  is  used,  factor  should  be  input  as  1.15 

28 

Creep  factor 

Multiplier  for  creep  curve,  should  be  based  on  bounds  being 
applied  to  creep  as  specified  in  ETL  1 1 10-2-365.  If  15% 
decrease  is  used,  factor  should  be  input  as  0.85 

29 

Initial  shrink,  strain 

Input  as  0.0 

30 

Reference  time 

Identifies  when  the  material  will  begin  gaining  strength  in  the 
analysis  (e.g.  use  5.0  if  material  being  described  is  in  the 

2nd  lift  and  5.0  day  placement  intervals  are  being  used) 

31 

Agg.  size  parameter 

Input  as  0.0 

32 

Reinforcement  ratio 

Input  as  0.0 
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6  Lock  and  Dam  Constitutive 
Modei  Input 


The  use  of  ABAQUS  or  ANACAP-U  subroutines  is  straightforward.  A 
portion  of  an  input  file  for  a  NISA  of  the  Olmsted  Lock  is  given  below  to 
clarify  the  required  ABAQUS  information  and  commands  when  using  a  self- 
defined  constitutive  model.  The  mathematically  derived  constants  that  were 
presented  in  the  chapters  regarding  creep,  autogenous  shrinkage,  and  aging 
modulus  are  contained  within  the  constitutive  model  subroutine  and  are  not 
input  variables.  The  variables  described  within  Table  1  are  typically  the  vari¬ 
ables  required  within  an  input  file.  The  ABAQUS  commands  and  Olmsted 
values  for  these  values  are: 
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*SOLID  SECTION, 

MATERIAL=M  1  ,ELSET=LIF1 
*MATERIAL,NAME=M  1 
*USER  MATERIAL, 
CONSTANTS=32 

4.675.. 100.E-6,2.1E6,.15,4.E-6,66.6,1 
0.,0.,0.,0.,0.,0.,0.,0. 
0.,0.,0.,0.,0.,0.,0.,0. 

2,0.4999,1. 1,0.9,0.,0.,0.,0. 

*DEPVAR 

23 

*SOLID  SECTION, 
MATERIAL=M2,ELSET=LIF2 
*MATERIAL,NAME=M2 
*USER 

MATERIAL,CONSTANTS=32 

4.675.. 1.)0.E-6,2.1E6,.15,4.E-6,66.6,1 
0.,0.,0.,0.,0.,0.,0.,0. 
0.,0.,0.,0.,0.,0.,0.,0. 

2,0.4999,1. 1,0.9,0.,10.,0.,0. 

*DEPVAR 

23 

*SOLID  SECTION, 
MATERIAL=M3,ELSET=LIF3 
*MATERIAL,NAME=M3 
♦USER  MATERIAL, 
CONSTANTS=32 

4.675.. 100.E-6,2.1E6,.15,4.E-6,66.6,1 
0.,0.,0.,0.,0.,0.,0.,0. 
0.,0.,0.,0.,0.,0.,0.,0. 

2,0.4999,1. 1,0.9,0.,20.,0.,0. 

♦DEPVAR 

23 

♦SOLID  SECTION, 
MATERIAL=M4£LSET=LIF4 
♦MATERIAL,NAME=M4 
♦USER  MATERIAL, 
CONSTANTS=32 

4,675.,100.E-6,2.1E6,.15,4.E-6,66.6,1 

0.,0.,0.,0.,0.,0.,0.,0. 

0.,0.,0.,0.,0.,0.,0.,0. 

2,0.4999,1. 1,0.9,0.,30.,0.,0. 

♦DEPVAR 

23 


Definition  for  lift  No.  1,  element  set 
LIFl 

Identifies  material  as  a  user-defined 
material 

Input  parameters  needed  for  user- 
defined  model  as  described  in 
Table  1 

Identifies  state-dependent  variables 

Definition  for  lift  No.  2,  element  set 
LIF2 


Definition  for  lift  No.  3,  element  set 
LIF3 


Definition  for  lift  No.  4,  element  set 
LIF4 
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7  Cracking  Criteria  and 
Model 


A  major  reason  for  performing  a  NISA  of  a  mass  concrete  stmcture  is  to 
predict  the  potential  for  cracking  to  occur  during  the  construction  and  service 
life  of  the  stmcture.  Finding  the  crack  potential  within  a  massive  concrete 
stmcture  provides  the  designer  with  important  information  regarding  the  qual¬ 
ity  of  the  stmcture.  If  the  analysis  indicates  a  high  potential  for  cracking 
during  constmction,  the  constmction  procedures  could  be  modified,  the  mate¬ 
rial  constiments  could  be  changed,  or  the  stmctural  geometry  could  be  rede¬ 
signed  in  order  to  reduce  this  potential,  thereby  providing  a  better  product  and 
a  more  reliable  stmcmre.  In  order  to  check  this  potential,  an  accurate  repre¬ 
sentation  of  the  principal  strains  are  needed  coupled  with  an  accurate  cracking 
criterion.  ANACAP-U  is  a  subroutine  that  provides  an  age-dependent  cracking 
criterion  to  be  used  in  conjunction  with  the  constitutive  model. 

Regions  of  potential  cracking  are  stmcture,  climate,  material,  and 
constmction-procedure  dependent,  but  several  typical  situations  can  be 
addressed.  Stmctural-related-cracking  potential  is  generally  higher  at  comers 
and  abmpt  changes  in  geometry.  Qimate-related-cracking  potential  is  gener¬ 
ally  higher  for  periods  of  cold  weather  during  and  after  constmction.  Material- 
related-cracking  potential  is  heavily  dependent  on  the  concrete  constituents 
affecting  the  adiabatic  temperature  rise,  the  aging  modulus,  the  creep  proper¬ 
ties,  and/or  the  shrinkage  properties.  Constituents  that  can  have  a  pronounced 
effect  are  aggregate  size,  fly  ash  content,  cement  type,  and  mixtures. 
Constmction-related-cracking  potential  is  heavily  dependent  on  procedures 
such  as  placement  ternperamre,  insulation,  lift  height,  lift  placement  rate,  lift 
placement  seasons,  and/or  lift  placement  sequence. 

The  potential  for  cracking  at  any  integration  point  in  a  finite  element  grid  is 
checked  using  an  interactive  stress-strain  cracking  criterion.  The  cracking 
criterion  is  not  explicitly  time  dependent,  which  is  why  an  interactive  stress- 
strain  criterion  is  used  where  the  time  effects  are  accounted  for  through  the 
age-dependent  modulus.  If  the  cracking  criterion  is  violated,  a  crack  will  be 
introduced  perpendicular  to  the  direction  of  the  maximum  principal  strain.  If  a 
crack  is  introduced,  the  constimtive  matrix  is  reformulated  within  the 
ANACAP-U  constimtive  subroutine  and  a  new  stress  state  is  developed  based 
on  zero  stress  in  the  principal  tensile  strain  direction.  The  new  constimtive 
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matrix  and  stresses  are  then  used  for  subsequent  calculations  until  another 
crack  is  indicated  by  the  criterion  or  until  the  crack  closes.  The  cracks  can 
close  when  placed  in  a  compressive  state  and  the  material  will  at  that  time  be 
capable  of  carrying  compressive  loads.  Depending  on  the  severity  of  the 
crack,  the  shear  resistance  is  reduced  at  the  cracked  integration  points,  but  the 
crack  wiU  have  limited  shear  resistance  which  is  a  function  of  friction  and 
aggregate  interlock. 

The  cracking  criterion  is  both  stress  and  strain  dependent.  Figure  6,  a 
general  plot  of  experimental  fracture  data  for  concrete,  shows  that  fracture  can 
occur  in  instances  of  high  stress-low  strain  or  high  strain-low  stress.  The 
triaxial  tension  and  cylinder  split  tests  are  of  the  high  stress-low  strain  variety 
while  the  uniaxial  compression  split  test  is  of  the  high  strain-low  stress  variety. 
The  latter  case  is  easily  explained  using  a  cube  loaded  in  compression  on 
opposite  faces,  principal  stress  direction  2  (assuming  that  principal  stresses  are 
arranged  from  highest  tensile  to  highest  compressive).  Ultimately,  the  cube 
wiU  fracture  due  to  the  strain  in  principal  stress  direction  1.  The  stress  in  this 
direction  will  be  small  tensile  stress,  which  is  an  indication  that  fracture  is 
related  to  the  strain.  The  strain  2^  would  be  -vej  which  would  imply  that  if 
the  compressive  stress  is  approaching  the  ultimate  compressive  stress,  the 
cracking  strain  would  be  v  times  the  uniaxial  ultimate  compressive  strain. 

Since  v  is  approximately  0.15  to  0.25,  the  fracture  strain  would  be  approxi¬ 
mately  15  to  25  percent  of  the  ultimate  compressive  strain.  (A  commonly 
assumed  value  is  10  percent  of  the  ultimate  compressive  strain.)  Therefore,  a 
strain-dependent  cracking  criterion  seems  reasonable.  The  first  case  is  self- 
explanatory.  Placing  a  cube  in  a  triaxial  state  of  tension  will  cause  failure 
prior  to  reaching  the  uniaxial  tension  failure  strain;  therefore,  it  is  also  highly 
stress  dependent. 

These  two  cases  indicate  that  a  criterion  which  is  dependent  on  stress  and 
strain  is  required.  In  addition,  the  time  dependency  is  needed  due  to  aging  of 
the  material  and  creep.  Therefore,  Figure  6  is  approximated  by  using  the 
fracture  strain  (obtained  from  a  slow  load  test),  the  fracture  stress 
(obtained  from  a  slow  load  test),  and  the  aging  modulus  E{%)  to  predict  the 
appropriate  stress-strain  failure  surface  defined  by  Of  and  Zj-  as  shown  in  Fig¬ 
ure  7.  The  slow  load  test  is  performed  on  a  simply  supported  beam  with  two 
point  loads  applied  in  a  manner  to  achieve  a  constant  moment  in  the  center  of 
the  beam.  The  load  is  increased  weekly  to  achieve  an  increase  of  tensile  stress 
in  the  extreme  fibers  of  25  psi.  This  procedure  is  repeated  until  the  beam 
cracks,  providing  the  slow  load  fracture  strain  e^,  the  age  of  the  concrete  x,  and 
the  slow  load  fracture  stress  a^.  Previous  tests  on  concrete  cylinders  using  the 
same  mixture  as  the  beam  provide  the  aging  modulus  £(x),  which  is  the  last 
piece  of  information  required  to, generate  the  failure  surface  in  Figure  7. 

The  cracking  criterion  is  either:  yes,  the  material  has  cracked,  or,  no,  it  has 
not.  This  yes/no  crack  prediction  is  necessary  and  correct  when  finding  the 
ultimate  response  of  the  structure,  but  it  is  not  very  useful  in  predicting  reli¬ 
ability  or  potential  for  cracking.  Therefore,  the  subroutine  developed  by  the 
ANATECH  Research  Corp.  provides  a  percentage  of  the  cracking  criterion  in 
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Figure  6.  Approximation  of  stress-strain  failure  curve  based  on  laboratory  test 


order  to  evaluate  the  potential  for  cracking.  A  percentage  approaching  100 
indicates  an  increasing  possibility  of  cracking.  Any  structure  with  a  NISA  that 
indicates  a  high  potential  for  cracking  as  described  in  ETL  1110-2-365  should 
be  evaluated  for  the  severity  of  the  consequences  of  the  predicted  cracks.  If 
the  consequences  are  deemed  detrimental  with  respect  to  either  safety  or  eco¬ 
nomics,  the  structure  should  be  redesigned  to  mitigate  the  effects  or  potential 
of  the  cracking.  Possibilities  for  redesign  include  but  are  not  limited  to  the  use 
of  additional  reinforcement,  the  revision  of  construction  procedures,  and/or  the 
modification  of  the  material  constituents  to  alleviate  or  control  the  cracking. 


The  cracking  criterion  is  a  stress-strain  interactive  criterion.  It  accounts  for 
age  dependency  of  the  criterion  through  the  use  of  a  linear  relationship 
between  cracking  stress  and  cracking  strain  which  is  dependent  on  the  aging 
modulus  as  shown  in  Figure  8.  With  the  appropriate  values  for  stress  and 
strain  at  a  given  time.  Figure  8  can  be  used  to  decide  if  a  region  of  the  struc¬ 
ture  has  cracked.  If  the  principal  stresses  and  their  respective  principal  strains, 
when  plotted  in  Figure  8,  are  within  the  triangle  enclosed  by  the  failure  surface 
and  the  two  axes,  no  cracking  occurs,  and  the  cracking  potential  is  calculated. 
If  the  point  of  principal  stress  versus  principal  strain  lies  outside  the  triangle, 
the  concrete  has  cracked.  If  the  system  is  assumed  to  be  cracked,  the 
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Stress 


Figure  7.  Cracking  failure  surface  and  e,  generation  from  the  slow  load  fracture  test  data 


Figure  8.  Computer-generated,  time-dependent  (aging  modulus)  cracking  failure  surfaces 
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constitutive  matrix,  stress  state,  nodal  forces,  and  stiffness  matrix  are  adjusted 
prior  to  continuation  of  the  analysis. 

The  failure  surface  is  a  function  of  the  fracture  stress  from  the  slow  load 
test  0^,  the  fracture  strain  from  the  slow  load  test  e^,  and  the  aging  modulus  at 

the  time  of  fracture  £/t),  where  x  is  the  age  of  the  concrete.  The  strain  axis 
intercept  is  determined  as 


£/  =  E.  ^ 


Eix) 


(23) 


as  shown  in  Figure  7.  This  intercept  value  remains  constant  for  the  entire 
NISA  and  is  a  prediction  of  the  concrete  cracking  strain  at  zero  stress. 
ABAQUS  input  data  require  the  user  to  input  a  cracking  strain  of 


(24) 


All  of  these  data  are  easily  obtained  from  a  slow  load  fracture  test.  The 
factor  of  1/2  is  a  function  of  the  input  required  by  the  ANATECH-developed 
subroutine  used  to  generate  the  correct  strain  axis  intercept  used  for  checking 
the  cracking  criterion.  Since  the  strain  intercept  remains  constant,  the  age 
dependency  is  related  to  the  time  variation  of  the  aging  modulus  (Figure  8)  for 
the  region  of  the  structure  being  evaluated  for  cracking.  This  concept  is  illus¬ 
trated  in  Figure  8.  The  stress  axis  intercept  for  a  given  age  x,  is  determined  as 


(25) 


Figure  8  shows  three  different  failure  surfaces  for  the  concrete  ages  of  Xj,  Xj 
and  Xj. 

Cracking  potential  is  a  quantitative  measure  of  the  imminence  of  violating 
the  cracking  criteria.  It  is  equivalent  to  the  ratio  of  1,  to  the  total  length  (Ij  -t- 
Ij),  as  shown  in  Figure  9,  where  1,  is  the  distance  from  the  origin  to  the  point 
(e,  c)  which  reflects  the  actual  principal  stress  and  strain  in  the  region  of  the 
structure  being  considered  for  cracking.  The  value  (1,  -t- 12)  is  the  length  of  the 
line  from  the  origin  to  the  failure  surface  passing  through  (e,  a),  which  are  the 
strain  and  stress  calculated  by  the  program.  The  cracking  potential  is  an  indi¬ 
cator  of  how  near  the  current  stress-strain  state  for  a  given  integration  point  is 
to  the  cracking  surface. 

The  following  is  the  algorithm  used  by  the  UMAT  subroutine  to  check  for 
cracking  as  well  as  the  process  which  occurs  when  cracking  does  and  does  not 
occur. 
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Figure  9.  Cracking  potential  generation  for  a  specific  cracking  failure  surface 


a.  Plot  the  point  represented  by  the  maximum  tensile  principal  stress  Cj 
and  its  respective  principal  strain  Ej  ;  check  if  the  point  is  inside  or 
outside  the  failure  surface. 

b.  If  inside  the  surface,  no  cracking  occurs,  the  cracking  potential  is  cal¬ 
culated,  and  the  next  integration  point  is  checked. 

c.  If  on  or  outside  the  surface,  introduce  a  crack  perpendicular  to  the 
direction  of  the  maximum  tensile  principal  strain. 

d.  Then  the  stress  in  this  direction  must  be  set  to  zero,  and  the  other  stres¬ 
ses  must  be  modified  to  reflect  that  change. 

e.  The  stiffness  matrix  must  then  be  modified  to  reflect  zero  load-carrying 
capabilities  in  that  direction  until  the  crack  closes  and  enters  a  compres¬ 
sive  state. 

/.  If  the  material  enters  a  compressive  state,  the  crack  is  assumed  to  have 
closed  and  100  percent  of  the  compressive  stiffness  is  reinstated  in  the 
direction  perpendicular  to  the  crack.  Once  the  material  is  placed  in  a 
tensile  state  again  the  crack  and  a  zero  stress  state  is  reintroduced  at 
this  location. 


Chapter  7  Cracking  Criteria  and  Model 


27 


The  3-D  case  requires  the  use  of  three  components  of  principal  stress  and 
strain  to  determine  the  cracking  potential  within  the  massive  concrete  structure. 
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8  Summary 


The  use  of  an  accurate  constitutive  model  is  imperative  when  performing  a 
mass  concrete  NISA.  ANATECH  Research  Corporation’s  ANACAP-U,  consti¬ 
tutive  relationship,  software  coupled  with  the  general-purpose  finite  element 
program  ABAQUS  provide  the  means  for  performing  reliable  NlSAs.  The 
ANACAP-U  constitutive  model  equations  can  be  adapted  to  fit  the  experimen¬ 
tal  data  for  a  given  project,  and  this  experimentally  based  constitutive  model 
can  then  be  implemented  with  ABAQUS’  user-defined  material  property  sub¬ 
routines  to  accurately  analyze  the  mass  concrete  structure.  The  creep,  shrink¬ 
age,  aging  modulus,  and  cracking  criteria  are  modeled  within  ANACAP-U  and 
are  considered  the  minimum  necessary  properties  for  accurate  modeling  of  the 
material  behavior.  The  general  theories  used  within  the  constitutive  model 
have  been  presented  within  this  report  in  a  practical  manner  to  enhance  the 
user’s  knowledge  regarding  the  use  of  material  constitutive  law  modeling  and 
its  use  in  analyzing  structural  behavior.  This  practical  discussion  is  far  from 
being  all-inclusive  and  users  of  the  ANACAP-U  constitutive  model  should 
read  its  manual  carefully  and  only  deviate  from  that  manual  after  obtaining  a 
thorough  understanding  of  the  software  and  material  modeling  concepts.  The 
use  of  these  constitutive  models  and  the  development  of  the  appropriate  con¬ 
stants  within  the  equations  representing  the  material  properties  should  be  a 
joint  or  team  effort  including  the  materials,  structural,  and  construction  engi¬ 
neering  staff  for  the  project.  The  intent  of  performing  a  NISA  is  to  provide 
more  reliable  and  cost-effective  structures  and  the  use  of  an  accurate  constitu¬ 
tive  or  material  model  is  necessary  for  this  intent  to  become  a  reality. 
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Appendix  A 

Calibration  of  Creep  Curve  to 
Data — Example 


Calibration  of  the  creep  curve  is  somewhat  more  involved  than  calibration 
of  the  modulus  of  elasticity  and  the  shrinkage  curves  since  the  equations  in  the 
constitutive  model  for  creep  are  not  directly  correlated  to  the  test  data.  The 
modulus  and  the  shrinkage  can  be  computed  with  the  equations  and  compared 
directly  with  the  test  results.  For  creep,  once  an  equation  has  been  defined,  a 
numerical  simulation  of  the  test  must  be  performed.  The  numerically  gener¬ 
ated  results  are  then  compared  with  the  plots  of  the  test  results.  As  described 
in  the  main  body  of  this  report,  steps  can  be  taken  to  obtain  a  first  approxima¬ 
tion  of  the  equation  of  the  creep  curve,  but  this  will  rarely  provide  the  final 
parameters  needed  to  satisfy  a  reasonable  calibration  process. 

In  order  to  obtain  an  accurate  equation  for  creep,  an  iterative  process  is 
usually  required  which  involves  changing  a  parameter  in  the  equation,  per¬ 
forming  the  numerical  analysis  of  the  test,  and  comparing  the  results  with  the 
test  data.  After  viewing  the  results,  this  process  may  need  to  be  performed 
again.  There  is  no  specific  number  of  iterations  required  and  they  will  vary 
from  one  concrete  mixture  to  the  next.  This  appendix  was  assembled  to  dem¬ 
onstrate  the  process  that  is  required  to  obtain  a  fit  of  the  creep  curve  equation 
to  the  test  data. 

The  test  data  used  in  the  calibration  process  presented  in  this  appendix  are 
for  the  interior  mixture  from  the  McAlpine  Lock  Replacement  project.  A  full 
set  of  tests  was  conducted  for  the  McAlpine  project  as  part  of  the  material 
characterization  study  (1994,  DM  No.  1,  Concrete  Materials,  McAlpine  Lock 
Replacement)'  and  for  the  nonlinear,  incremental  structural  analysis  (NISA) 
that  was  performed  (in  preparation,  Fehl,  Riveros,  and  Gamer). 

Using  the  process  to  define  the  creep  equation  as  outlined  in  the  main  body 
of  the  report,  3-day  creep  data  were  selected  as  the  basis  for  establishing  the 
parameters.  The  times  selected  for  curve  fitting  were  4,  20,  and  100  days. 


^  References  cited  in  this  appendix  are  located  at  the  end  of  the  main  text. 
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Using  these  times  and  assuming  the  exponential  terms  of  Equation  4  of  the 
main  text  are  equal  to  a  value  of  0.001,  values  for  Tj,  r2,  and  can  be  com¬ 
puted  as  follows; 


e  =  0.001  ^  r, 


^  -^3(100-3) 


0.001  => 


6.9078 

0.4063 

0.0721 


For  this  case  the  constant  term  in  Equation  4  was  replaced  with  a  third 
exponential. 

Now  that  values  for  rj,  r2,  and  have  been  computed  they  can  be  substi¬ 
tuted  into  the  equation.  The  values  for  J(t,x)  for  times  r  of  4,  20,  and  100  days 
and  T  of  3  days  can  be  found  from  the  test  data  in  the  material  characterization 
study  report  (USAED  Louisville  1993)  and  substimted  into  the  equation. 

Values  from  the  Me  Alpine  project  study  were: 

7(4,3)  =  0.2  millionths/psi 

7(20,3)  =  0.6  millionths/psi 

7(100,3)  =  0.8  millionths/psi 

So,  once  the  values  for  r^,  r2,  and  and  7(/,t)  are  substituted,  then  three  equa¬ 
tions  with  three  unknowns  remain  and  can  be  solved  simultaneously.  The 
three  equations  are: 


0.2  X  10  ®  =  0.999Ai 

-1-  O.3339A2 

+  O.O687A3 

(Al) 

0.6  X  10  ®  =  Aj 

+  O.999A2 

+  O.298OA3 

(A2) 

0.8  X  10  ®  =  Aj 

+  A2 

+  O.999A3 

(A3) 

Solving  the  equations  simultaneously  results  in  the  following  values  for  Aj,  A2, 
and  A3: 

Ai  =  0.12513  X  10  ® 

A2  =  0.3904  X  10  ® 

A3  =  0.28475  X  10-® 

These  values  along  with  those  computed  above  for  r,,  r2,  and  r-^  can  be  used  in 
the  equation  for  creep  in  the  constitutive  model  and  a  numerical  analysis  of  the 
creep  cylinder  tests  can  be  performed.  The  strains  from  these  analyses  can  be 
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be  compared  with  the  actual  strains  recorded  during  the  test  to  determine  if  the 
parameters  computed  are  adequate. 

The  analyses  were  performed  using  these  initial  parameters  and  the  results 
from  the  analyses  are  compared  with  the  test  results  in  Figure  A1  for  the 
1-day,  3-day,  and  14-day  creep  tests.  The  magnitudes  are  approximately  cor¬ 
rect  for  the  3-  and  14-day  tests  but  the  shape  of  the  curves  is  not  satisfactory 
for  these  two  tests.  Adjustments  must  be  made  and  then  the  analyses  per¬ 
formed  again  and  the  results  compared  with  the  test  results. 

Magnitudes  of  the  curves  are  generally  controlled  by  the  factors  while  the 
shapes  of  the  curves  can  be  changed  by  adjusting  the  exponentials.  Since  the 
results  from  the  analysis  are  decreasing  compared  with  the  test  results,  the 
exponential  for  the  long  term  (r^)  will  be  adjusted  as  shown  in  Table  A1  for 
attempt  No.  2.  Results  of  the  analysis  using  values  for  attempt  No.  2  and 
results  from  attempt  No.  1  and  from  the  test  data  are  shown  in  Figure  A2. 

The  long-term  shapes  are  improved  but  the  results  are  too  low  for  the  3-  and 
14-day  tests. 

Since  the  early  time  drop  appears  to  be  excessive,  the  exponent  for  early 
times  will  be  reduced  to  change  the  behavior  of  the  early  time  curve.  Attempt 
No.  3  will  revise  exponent  as  shown  in  Table  Al.  The  results  of  attempt 
No.  3  are  shown  in  Figure  A3  and  are  compared  with  the  results  from  attempt 
No.  2  and  the  test  data.  The  revised  exponent  changed  the  early  time  behavior 
of  the  curve  slightly  and  increased  the  strain  recorded  for  aU  three  tests,  but 
the  values  are  still  low  and  the  increase  in  strain  is  still  too  gradual. 

To  increase  the  steepness  of  the  curve  in  the  first  50  days,  the  exponential 
term  will  be  reduced  further  as  shown  in  Table  Al.  Exponent  r-^  will  also 
be  changed  to  flatten  the  long-term  curve.  The  results  of  attempt  No.  4  are 
shown  in  Figure  A4  and  are  compared  with  the  results  of  attempt  No.  3  and 
the  test  data.  In  all  three  tests  the  curve  is  now  approaching  a  shape  compara¬ 
ble  to  that  of  the  test  data,  but  the  magnitude  of  the  strains  is  too  low  for  the 
3-  and  14-day  tests. 

To  increase  the  magnitude  of  the  strains,  the  A2  term  will  be  increased  to 
revise  the  overall  magnitude  of  the  curve.  The  change  made  for  attempt  No.  5 
can  be  seen  in  Table  A 1,  and  the  results  of  the  analyses  are  shown  in  Fig¬ 
ure  A5  and  are  compared  with  results  from  attempt  No.  4  and  the  test  data. 

As  can  be  seen  in  Figure  A5  the  numerical  results  now  exceed  the  values  from 
the  test  data  for  all  three  tests. 

Two  parameters  will  be  changed  on  the  next  attempt.  The  early  time  factor 
Aj  win  be  increased  substantiaUy  to  maintain  the  steep  curve  at  early  times 
and  factor  A3  wiU  also  be  decreased  significantly  to  try  and  maintain  a  flat 
curve  at  later  times.  The  changes  made  are  shown  in  Table  Al  and  the  results 
shown  in  Figure  A6.  As  can  be  seen  in  Figure  A6  the  change  in  parameters 
basicaUy  decreased  the  steepness  of  the  slope  and  the  magnitudes  are  stiU  too 
high  at  longer  times. 
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Table  A1 

Exponentials  and  Factors  Used  In  Calibrating  Creep  Curve 


Exponentials 


Factors 


Attempt  No. 


1 

6.9078 
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6.9078 
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t.O 

4 

O.t 

5 

0.1 
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0.1 
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1.0 
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1.0 
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1.0 

10 

1.0 

0.4063 


0.4063 


0.4063 


0.4063 


0.4063 


0.4063 
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002 
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003 
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Figure  A2.  Time  history  strain  plots  of  creep  test  data  and  numerical  results  for 
attempts  No.  1  and  2  (Continued) 
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c.  14-day  test 


Figure  A2.  (Concluded) 


Since  steepness  of  the  curves  again  appears  to  be  a  problem,  the  exponen¬ 
tial  terms  will  again  be  adjusted.  This  time  all  three  exponentials  will  change 
as  shown  in  Table  A1  with  resulting  plots  as  shown  in  Figxire  A7.  The 
changes  made  increased  the  early-time  steepness  of  all  three  curves  signifi¬ 
cantly  but  did  result  in  a  reduction  in  magnitude  of  strains  at  later  times. 

Since  an  earlier  increase  in  (attempt  No.  6)  caused  a  significant  increase 
in  the  steepness  of  the  early-time  curve,  Aj  will  now  be  reduced  to  near  its 
original  value  as  shown  in  Table  Al.  To  account  for  the  loss  in  magmtude, 
factor  A3  will  be  increased.  Figure  A8  shows  the  results  and  as  can  be  seen 
the  curve  is  now  much  flatter  but  the  magnitude  is  too  low  except  for  the 
1-day  test.  In  addition,  curves  for  3  and  14  days  should  be  a  little  steeper  at 
early  times. 

To  get  the  order  of  magnitude  correct  for  the  3-  and  14-day  curves,  the 
factor  A2  will  be  increased  as  shown  in  Table  Al.  The  effect  of  this  increased 
factor  is  shown  in  Figure  A9  where  curves  for  3  and  14  days  are  getting  very 
close  to  the  test  values.  Magnitudes  are  still  a  little  high  on  the  3-day  curve 
and  the  steepness  of  the  curve  over  the  long  term  should  be  decreased  slightly. 
Values  for  the  1-day  creep  are  too  high,  but  it  is  becoming  apparent  that  expo¬ 
nents  and  factors  that  can  be  used  to  fit  the  3-  and  14-day  curves  will  not  pro¬ 
vide  good  agreement  with  the  1-day  curve. 
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Figure  A3.  (Concluded) 


The  tenth  attempt  will  increase  the  exponent  to  flatten  the  curve  some¬ 
what  and  the  factor  A3  will  be  decreased  to  reduce  the  magnitude  over  the  long 
term.  The  changes  are  shown  in  Table  A1  and  the  results  are  plotted  in  Figure 
AlO.  As  can  be  seen  in  Figure  AlO  the  magnitude  of  strains  decreased  for  aU 
three  tests  for  the  portion  of  the  curves  beyond  50  days.  Figure  A  10(b)  shows 
that  the  data  for  attempt  No.  10  provide  very  good  agreement  for  the  3-day 
data  and  reasonable  agreement  with  the  14-day  data  as  shown  in  Fig¬ 
ure  A  10(c).  As  stated  in  the  paragraph  above,  a  discrepancy  exists  for  the 
1-day  creep  test  results  but  this  does  not  appear  to  be  able  to  be  resolved  with¬ 
out  ruining  the  fits  of  the  3-  and  14-day  curves. 

The  curve  fits  shown  for  attempt  No.  10  in  Figure  AlO  could  be  adjusted 
further  to  try  to  account  for  the  high  values  for  the  1-day  test  and  to  increase 
values  at  the  later  times  for  the  14-day  test.  A  decision  could  also  be  made  to 
select  attempt  No.  9  factors  and  exponentials  since  these  values  provide  an 
excellent  fit  with  the  14-day  creep  data.  The  difficulty  always  with  calibrating 
the  creep  curve  though  is  that  a  change  made  to  correct  a  deficiency  in  the  fit 
for  one  curve  will  affect  the  other  two  curves.  Although  difficult,  factors  and 
exponentials  can  be  found  that  will  accurately  reflect  aU  three  sets  of  the  test, 
but  generally  fitting  the  data  to  two  of  the  curves  will  have  to  suffice. 

As  can  be  seen  by  the  steps  taken  in  this  appendix  to  calibrate  a  creep 
curve,  the  process  can  be  tedious.  The  process  was  presented  to  provide  the 
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Figure  A4.  Time  history  strain  plots  of  creep  test  data  and  numerical  results  for 
attempts  No.  3  and  4  (Continued) 
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c.  14-day  test 


Figure  A4.  (Concluded) 


reader  a  better  understanding  of  how  involved  the  process  of  calibrating  the 
creep  equation  can  be  and  to  provide  insight  on  how  to  adjust  the  parameters. 
Calibration  of  the  creep  curve  is  certainly  not  exact  but  through  careful  manip¬ 
ulation  of  the  parameters  and  use  of  good  engineering  judgment  a  set  of 
parameters  can  be  found  which  adequately  capture  the  creep  behavior. 
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Figure  A5.  (Concluded) 
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Figure  A6.  Time  history  strain  plots  of  creep  test  data  and  numerical  results  for 
attempts  No.  5  and  6  (Continued) 
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c.  14-day  test 


Figure  A7.  (Concluded) 
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Figure  A9.  Time  history  strain  plots  of  creep  test  data  and  numerical  results  for 
attempts  No.  8  and  9  (Continued) 
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Figure  A10.  (Concluded) 
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